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We present a calculation of the two-photon and two-gluon widths for the equal mass quarkonium 
states 1 So, 3 Po and 3 P2 of the charmonium and upsilon systems. The approach taken is based 
on using the full relativistic qq — > 77 amplitude together with a wave function derived from the 
instantaneous Bethe-Salpeter equation. Momentum space radial wave functions obtained from an 
earlier fit of the charmonium and upsilon spectra are used to evaluate the necessary integrals. 
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1. INTRODUCTION 

Equal-mass quarkonia are eigenstates of the charge conjugation operator C with eigenvalues C = (— 1) L+S . 
As such, the 1 Sq, 3 Pq and 3 P2 levels of charmonium and the upsilon system can decay into two photons. 
These same states can also decay into two gluons, which accounts for a substantial portion of the hadronic 
decays for states below the cc or bb threshold. One of us (WWR) was previously involved in an investigation 
but, as more data become available, there has been a renewed interest in this subject 



In the calculation of the two gamma widths that follows, we wish to include the relativistic and QCD 
effects in the wave functions that are used to compute the two photon decay amplitudes. We do this by 
making use of the variational wave functions obtained in [9[. These wave functions were computed using 
a semi-relativistic model containing both the v 2 /c 2 and one-loop QCD corrections to the potential and 
optimized to provide an accurate description of the cc and bb quarkonium spectra. One of the parameters 
determined in this process is the renormalization scale appropriate to the particular quarkonium system. 
We use the resulting radial functions to construct the individual 1 Sq, 3 Pq and 3 P2 wave functions with the 
proper spin dependence obtained from a decomposition of the instantaneous Salpeter equation. One of the 
objectives of taking this approach is to investigate how the inclusion of the full wave function information 
compares with the practice of using the square of the radial wave function |i? n o(0)| 2 for s-statcs or |-R^i(0)| 2 
for p-states. This much can be accomplished by using the expression for the invariant amplitude, M, which 
can be derived using the instantaneous Salpeter wave function, 4>{p) 

M =e 2 Jd 3 P Tr [C^^Sp^-k)^^ + C~ 1 ^S F [p-h!)^*<i>{p)] . (1) 

Here, C is the charge conjugation matrix, k and k' are the photon momenta, e and e' are the photon 
polarization vectors, Sf(p) is the quark propagator and p^ — (p, i\Jp 2 + to 2 ) = (p,iE), with m denoting 
the quark mass. The wave function 4>(p) is the Fourier transform of the instantaneous position space wave 
function an d its radial portion is obtained from the variational calculation in Ref.[9(. 
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There are also QCD corrections to the two gamma widths. Because of color conservation, these corrections 
are finite at the one-loop level. They could, in principle, be calculated by extending Eg. ([I]) as a perturbativc 
series in the QCD coupling as, but the one- loop corrections are known EH 13 EH and will be included 
in the final results. 

In the next section, we construct the necessary wave functions, evaluate the trace and calculate the widths 
T 77 for the 1 So, 3 -Po and z Pi states of the cc and bb systems. We conclude with a numerical evaluation of 
the two-photon widths and some comments on the two-gluon widths T gg . 



2. EVALUATION OF THE BOUND-STATE DECAYS 



The wave functions for the solution of the instantaneous Salpeter equation can be decomposed into spin 
singlet and spin triplet forms, which are [l4| 
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for the singlet states and 
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for the triplet states. These states have the conventional normalization 



J d 3 P P*(p)P(p) = 1 J d 3 P V\p)-V(p) = 1 , 



where 



P(p) = 
V(p) = 

Vi(p) = 



-in 
1 



P<Pnl(p) 



for 



^f ] Mm(p) for 3 P 2 



The spin-two polarization vectors Qj satisfy Qj 



fjf 3 and df > 
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The explicit form of the qq — > 77 amplitude in Eq. ([I]) is 



A = e'ei 



-2p-k 



-2p-k' 



C. 



(2) 
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(4) 
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where e is the proton charge and e q the fractional quark charge. With the aid of the wave functions in Eq. ([5]), 
the traces in Eq. |T]) can be evaluated rather straightforwardly for the spin states. In the quarkonium rest 
frame, including the factors l/2w from the photon normalization and l/(27r) 3 ' 2 from the Fourier transform 
of the position space wave function, the results are, 



Mi 1 S ) = e 2 e 2 q 



1 m f 3 k-(e'* x e*) 



8tt 2 lu 
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[(p-jfc) 2 -E 2 \ 



(7) 
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for the 1 S'o state, and 



1 m 



1 [uj{p-k) 2 e'*-e* + 2Ep-e'*p-e*} 



pE 
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0nl(p) : 
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for the 3 Pq state. The evaluation of Ai( 3 P 2 ) is somewhat more tedious, yielding 



i 



(M) 



j »tt 2 lu 2 J 'p[(p-k) 2 -E 2 

[uj(p-k) 2 e'*-s* + 2Ee'*-pe*-p}piPj 



E(e'*pjE* -p + e'* -pe*pj) — uie'* -s* p-kkiPj 

nl(p) , 



E(E + m) 

where we have dropped a term involving PiPjPk that vanishes upon solid angle integration. 
To complete the calculation of the amplitude, we need to evaluate the angular integrals 

1 
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(10) 
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noting that angular integrals of this type with an odd number of p^s in the numerator vanish. The values of 
the coefficients Aq through B 3 are given in the Appendix. In terms of these coefficients, the T amplitudes, 
including a factor of V3 for color, are 



Mi 1 So) 
for the 0~ state and 

M( 3 P ) = 
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for the + state. As with the traces, the 3 P 2 integration over dQp is more complicated since it involves four 
Pi's in the numerator. The result is 



M( 3 Po) = e 2 e 2 — -^- 
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(15) 



The calculation of the widths can be performed using the formula for the decay of a particle with spin J 
into two photons, 



r„ ^ J-w 2 — - — V V [dQt\M\ 
77 16tt 2 2J+ 1 ^ ^ J kl 1 



(16) 



M=-J pol 
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We sum over the photon polarizations using the Coulomb gauge 

Y,e*&)<k) = 6 ij -kk J : ( 17 ) 

pol 

and over the spin-two projections using 

2 1 1 

E f'"^ = 2^^« + " 3^ m " ■ (18) 

M=-2 



The resulting two-gamma widths are 
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(19) 
(20) 
(21) 



where conservation of energy requires that the photon energy u> satisfies cu = M q q/2. The last factors in 
Eqs. Q19 p -(|2i p include the one-loop QCD correction. The two-loop QCD correction to the decay rate has 
been examined in [l5| and appears to be large. We have not attempted to include this correction in our 
evaluation of the decay amplitudes. 



3. RESULTS AND CONCLUSIONS 



The integrals that remain, Eqs. (|19p - (|2~lj) , were evaluated using both the perturbative and non-perturbative 
wave functions from [9j with the aid of Mathematica. Our results for the charmonium system are shown in 
Table U fl6l , [l7| and those for the upsilon system system are given in Table [Til 



TABLE I: The data denoted by * were taken from Ref. [l6| and those denoted by t were taken from Ref. [J_3j. The 
column labeled Tnr+qcd is the value of the width calculated using |i?„o(0)| 2 or |74i(0)| 2 along with the QCD 
correction. 



Decay 


Tpcrt 


Tnon — Pert 


Lnr+qcd 




Vc -> 77 
Xco -» 77 
Xc2 -> 77 

Vc -» 77 


5.09 keV 
2.02 keV 
0.46 keV 
2.63 keV 


3.48 keV 
2.12 keV 
0.19 keV 
1.96 keV 


13.1 keV 
5.35 keV 
1.55 keV 
10.5 keV 


7.2 ±0.9 keV* 
2.53 ± 0.45 keV f 
0.60 ± .08 keV f 
< 7.0 ± 3.5 keV* 



For charmonium there is reasonable agreement between our calculations and the experimental results, 
with some preference for the perturbative calculation. The results are much smaller than the non-relativistic 
widths including the one-loop QCD correction. The two-gamma decay rates for the upsilon system have not 
been measured and the r]b ground state has only recently been observed [l8[ . The trend of the widths obtained 
using the non-relativistic results modified by QCD corrections is again to be larger than our calculated widths. 
Note that, in this formalism, knowledge of the masses is not necessary for the calculation of their 77 
widths since the only mass that occurs in Eq. (fT9| is the quark mass, which we obtain from Ref. This 
fact could account for the close agreement of our r 77 ( 1 S , ) results with those of Ref.[|[, even though these 
authors use a very different potential to get their radial wave functions. 



5 



TABLE II: 



Decay 


Tport 


Triori — Pert 


Tnr+qcd 


Vb 


-* 77 


0.30 keV 


0.32 keV 


0.55 keV 


Xbo 


-> 77 


32.9 eV 


94.4 eV 


58.4 eV 


Xb2 


_> 77 


7.19 eV 


5.48 eV 


9.85 eV 


Vb 


_» 77 


0.14 keV 


0.15 keV 


0.20 keV 


Xbo 


-> 77 


34.1 eV 


94.5 eV 


68.3 eV 


Xb2 


-> 77 


7.59 eV 


5.57 eV 


11.5 eV 


Vb 


_> 77 


0.10 keV 


0.11 keV 


0.22 keV 



Finally, we can convert our results from two-photon decays to two-gluon decays by replacing 3a 2 by 
2a|/3 in Eqs. (|19p - ([21"]) . This decay process accounts for a substantial portion of hadronic decays for states 
below cc or bb threshold. There arc, however, significant radiative corrections as well as contributions from 
three-gluon decays and, thus, the two gluon mode does not tell the whole story. Our results for these two- 
gluon decays are given in Tables IIIII and IIVI From Table III, it can be seen that the two-gluon widths are 
smaller than the hadronic widths of the charmonium states. 



TABLE III: 



Decay 


r ' Pert 


Triori — Pert 


Texp 


Vc -> gg 
Xco -> gg 
Xc2 -> gg 

v'c -> gg 


15.70 MeV 
4.68 MeV 
1.72 MeV 
8.10 MeV 


10.57 MeV 
4.88 MeV 
0.69 MeV 
5.94 MeV 


26.7 ±3.0 MeV 
10.2 ±0.7 MeV 
2.03 ±0.12 MeV 
14 ± 7 MeV 



TABLE IV: 



Decay 




r non — Pert 


Vb 


-> gg 


11.49 MeV 


12.39 MeV 


Xbo 


-> gg 


0.96 MeV 


2.74 MeV 


Xb2 


-> gg 


0.33 MeV 


0.25 MeV 


Vb 


-> gg 


5.16 MeV 


5.61 MeV 


Xbo 


-> gg 


0.99 MeV 


2.74 MeV 


Xb2 


-> gg 


0.35 MeV 


0.26 MeV 


Vb 




3.80 MeV 


4.11 MeV 
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APPENDIX A: ANGULAR INTEGRAL COEFFICIENTS 

The coefficients Aq through B3 occurring in Eqs. (|10p - (| 1 2[) are presented below [20| . 
Ep \E+pJ ' 
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TT 
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TT 



(3E 2 -p 2 )ln 



i 2 In 
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E+p 
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E + p 

2Ep 
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TT 



(105S 4 -90£;V + 9p 4 )ln 
(15£ 4 - + 3p 4 )ln 



E-p 
E + p 
E-p 
E + p 



\2Ep 



3to 4 In ( — 




\ E 


+ Pj 



{6E 3 p - lQEp 3 ) 



(210E 3 p- UOEp 3 
- (30E 3 p - 2QEp 3 ) 
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